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models are the same as of General Relativity. So, current experiments are not applicable for such 
search therefore different methods of observation and higher accuracy are required. 
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I. INTRODUCTION 

A set of multidimensional gravity models beginning 
from the Kaluza-Klein one [J] result from the attempts 
to construct an unified field theory. As we live in a 
space-time with four noncompact dimensions any mul¬ 
tidimensional theory needs an appropriate effective four- 
dimensional limit consistent with the predictions of gen¬ 
eral relativity (GR) and the results of observations and 
experiments. 

String theory Q along with the loop quantum grav¬ 
ity I is a promising candidate for a quantum theory of 
gravity nowadays. Lovelock gravity ||J appeared to be a 
ghost-free four-dimensional low-energy effective limit of 
the string theory [|, 6:j: 

L = \/—g (R + q: 2 <S'“ + T ...) , (1) 

where S n is Euler characteristic of the n-th order. The 
leading and the most studied among them is the sec¬ 
ond order curvature correction given by a Gauss-Bonnet 
term S 2 = S GB = A,,/,/ IP' 1 ' 1 - 4 RijR ij + /{/ The ef¬ 
fective four-dimensional limit of the string theory also 
includes the scalar field that is the projection of the 
g io io component of the ten-dimensional string metric 
to the four-dimensional manifold. The Gauss-Bonnet 
term coupled to the scalar (dilatonic) field [M3 de¬ 
scribes the influence of the compact extra dimensions on 
the four-dimensional space-time. Therefore the Gauss- 
Bonnet theory with the dilaton scalar field serves as the 
effective four-dimensional limit of the string theory. 

Unlike the string theory the Randall-Sundrum model 
allows the only extra dimension to be large and even 
infinite [H], [I]| ■ This model considers four-dimensional 
branes with tension embedded into a five-dimensional 
space-time (bulk) that is assumed to have an AdSs ge¬ 


ometry. All the matter and the three fundamental inter¬ 
actions are localized on this brane except gravity which 
is allowed to propagate into the bulk along the extra non¬ 
compact dimension. Thus the Randall-Sundrum model 
contains the description of the four-dimensional space- 
time from the very beginning and therefore does not 
need any special theory serving as its effective four- 
dimensional limit. The Randall-Sundrum I (RSI) model 
includes two branes with different properties helping to 
solve the hierarchy problem 0- Moving the second 
brane to the infinity led to the Randall-Sundrum II 
(RSII) model with one brane [lj|. In this paper we deal 
with RSII only. 

Black hole solution is a basic one for any theory of 
gravity. First of all, it describes the compact object that 
a very massive star at the end of its life cycle collapses 
into. It also features the curvature of the space-time 
produced by the presence of matter specific for the con¬ 
sidered gravity model. Any extended theory of gravity 
should be consistent with the predictions of GR and the 
observational results therefore the existence of black holes 
and their properties are important indicators of the the¬ 
ory’s viability. 

The Gauss-Bonnet solution has been studied explic¬ 
itly during the recent years @,§£3. On the contrary 
there are several different solutions for Randall-Sundrum 
model 0®. It was argued that static black holes 
cannot exist in RSII with a radius much greater than 
the AdS length l fl9l - [23 and that even very small RSII 
static black holes do not exist 0, [23} ■ Therefore RSII 
solutions for large black holes that have been found inde¬ 
pendently by Figueras and Wiseman [24], [2j| and Abdol- 
rahimi, Cattoen, Page and Yaghoobpour-Tari [26| are an 
important improvement of the Randall-Sundrum model 
interesting for further consideration. In their work [26} 
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Abdolrahimi, Cattoen et al. compare the obtained black- 
hole solution to the one Figueras and Wiseman [25] and 
show that these solution s ag ree closely. In this paper we 
use the ACPY solution [26| as it contains the necessary 
details. The Figueras-Wiseman solution 24] is consid¬ 
ered separately. 

The outline of this paper is as follows: in Sec¬ 
tion “Weak field limit” we discuss the weak field and 
slow motion approximation of the Gauss-Bonnet and 
Randall-Sundrum theories. Section “Thermodynamics 
and PBHs” is devoted to the analysis of thermodynami¬ 
cal properties of these models and their influence on the 
primordial black holes mass spectra. In Section “Discus¬ 
sion and Conclusions” we discuss the results obtained, 
offer conclusions, and outline the next steps. 


II. WEAK FIELD LIMIT 

As a weak-field limit we consider the dynamical con¬ 
ditions in the solar system i.e. post-Newtonian approxi¬ 
mation. The metric tensor g pv can be represented as the 
perturbation h pv around Minkowski space-time rj pu [27| : 


9^.v — V/iv Y (2) 

In this paper we consider only spherically-symmetric 
solutions and therefore the static gravity field at the 
distance r from its source. In the first post-Newtonian 
(PPN) order the correction to the gravitational field h 
can be expressed by the series with respect to the nega¬ 
tive powers of the radial coordinate r up to the following 
order terms: 

/loo ~ 0(r~ 3 ), hgj ~ 0(r~ 4 ), h tj ~ 0(r~ 2 ) (3) 

and use geometric units h = c = G = 1 with non- 
dimensional masses expressed in units of Plank mass. 

The PPN limit is well tested by experiments HU HU . 
The better the experimental accuracy becomes |28| the 
more opportunities to test small gravitational effects pre¬ 
dicted by currently viable theories should appear. We use 
the expansion d3j) to compare the magnitudes of the pre¬ 
dicted effects in order to see if the specific effects of the 
considered solutions can be tested. As the PPN requires 
the weak-field limit approximation we apply our results 
to solar system where the post-Newtonian parameters 
(PPN-parameters) are measured with high precision j32|. 
Our results are inapplicable to the strong field limit. 


A. Gauss-Bonnet gravity 

We begin our consideration by exploring the weak 
field limit of the Gauss-Bonnet theory (Here and further 
during solving the equations we use the dimensionless 
Planckian units, only at the step of numerical estimation 


values we jump to the usual ones) 




d A Xy/~^g 


R + 2d ll (j)d IJ '(l)+ 


+ Xe 2<t> ScB + ■ • ■ 


( 4 ) 


where (f> is the potential of the dilatonic field, Sgb is the 
Gauss-Bonnet term and A is the string coupling constant. 
For this purpose we construct a post-Newtonian param¬ 
eterization of the static asymptotically flat spherically- 
symmetric Gauss-Bonnet solution 

ds 2 = A dt 2 — dr 2 — r 2 (cLO 2 + sin 2 9d , 

A = l-^+0(r- 2 ), a = l + 0(r~ 2 ), 

4> = '/’oo + — + O (r 2 ), 

where (t, r, 9 , ip) are usual spherical coordinates and func¬ 
tions A and a depend on the radial coordinate r only, M 
is the Arnowitt-Deser-Misner (ADM) mass, D is dila¬ 
tonic charge i. e. the effective charge of the scalar field 
source and </>oo is the asymptotic value of dilatonic po¬ 
tential H[- As it was argued in [10] D oc 1/M. 

We put metric fl5J) with the expansions © to the field 
equations written down in the most computationally con¬ 
venient form |29j : 

= Stt (T™ + T+, + T £ B ) , ( 7 ) 

where T ™ is a matter stress-energy tensor while T// v and 
T?B describe the scalar field and the Gauss-Bonnet term 
presence: 


( 5 ) 

( 6 ) 


T /i = ^ d v <t> - i g pv d p </> dp^j , 


rp GB 

His 


^ [(V M V y -3^D)(e- 2 ^R) 

2 (□<£<£+ R P „) 

2V^V CT (e- 2 ^ R ppua ) . 


Using the standard computational techniques 27] the 
leading order for the nontrivial correction to the Gauss- 
Bonnet metric tensor is 


6h^=8— + 0(r~ 5 ). (8) 

Comparing this result with ([3| we see that the correc¬ 
tion term Q lies beyond the PPN order that should be 
proportional to 1/r 2 . Thus the parameters of the Gauss- 
Bonnet model cannot be constrained by the solar system 
tests. This result is consistent with conclusions from (2Sj 
where the cosmological limit of the discussed model was 
studied. 










3 


B. Randall-Sundrum gravity 


The black hole solution of the Randall-Sundrum model 
was constructed by Figueras and Wiseman [H[ using an 
associated five-dimensional anti-de Sitter space (AdSs) 
and dS 5 -CFT 4 correspondence (dfij . The Figueras- 
Wiseman solution describes a static black hole with ra¬ 
dius up to ~ 20£ and reproduces four-dimensional GR 
on the brane in the low curvature and the low energy 
limit. We intend to use the fact that Schwarchild metric 
can be used not only as a black hole one but also as a 
description of gravitational (stellar) system far from the 
central body. For example, for Solar system (taking into 
account all the limitations and corrections). 

The five-dimensional metric can be written near the 
conformal boundary 2 = 0 as 

l 2 

ds 2 = — \dz 2 + g^(z,x) dx^dx v \ , (9) 

where 2 is a coordinate of the brane along the extra di¬ 
mension and gfj, v (z,x) is the metric on the brane deter¬ 
mined by the Fefferman-Graham expansion [3C|. The 
corresponding effective four-dimensional field equations 
frit are: 


G 


Ills 


= 87T GO% ane + e 2 {l 6 ^G 4 (T M c / T [ 5 ]) + 
[g] +k>ge [s] | + 0(e 4 loge), 


( 10 ) 


where G 4 is the usual four-dimensional gravitational con¬ 
stant, T^l ane is the stress-energy tensor of the matter lo¬ 
calized on the brane, tensors {T^f T [ 5 ]), a^, v [g) and b^[g] 
result from the extra dimension and depend on the metric 
tensor components and e is a small perturbation param¬ 
eter indicating the deviation of the brane position from 
the equilibrium state z = 0 . 

The additional term in the post-Newtonian expansion 
of the Figueras-Wiseman solution calculated in this paper 
is 


Sh 


FW 

00 


121 e 2 M 2 


( 11 ) 


The obtained value m lies within the PPN limit © 
and points at a potentially observable effect. In Randall- 
Sundrum model gravity is allowed propagate into the 
bulk along the extra dimension therefore the effect de¬ 
scribed by HID most likely leads to the negative non¬ 
linearity in gravitational superposition. In other words 
the resulting gravitational field produced by two or more 
massive objects can be less than the direct vector sum of 
their contributions. The parameterized post-Newtonian 
(PPN) parameter /3 is responsible for such an effect 
(271 l28l |. Therefore the result m should be expressed 
as follows: 


where M is the mass of the massive central object. In 
the considered case it equals the solar mass. It also is 
expressed in Plank units of mass and therefore is dimen¬ 
sionless. 

The constraint on the PPN parameter /3 obtained from 
analysis of the lunar laser ranging data IH ^ m - 1 < 
1.1 x 10 4 [H]. The admitted region of the AdS length is 
limited by the results of the Newton’s law test i < 10 -5 m 
[33 |. Therefore the upper limit on the value of e is: 

e < 5.7 x 10" 47 cm « l Ph (13) 


Originally the parameter e was assumed to be negli¬ 
gibly small and the vanishing value found in (1131) im¬ 
plies that in fact e = 0. Thus the Figueras-Wiseman 
four-dimensional black hole solution is not only is self- 
consistent but well consistent with the solar system con¬ 
straints as well. Therefore this solution is indistinguish¬ 
able form GR in the PPN limit after all. 

The other recent Randall-Sundrum solution obtained 
by Abdolrahimi, Cattoen, Page and Yaghoobpour- 
Tari (ACPY) (26| is asymptotically conformal to the 
Schwarzschild metric and includes a negative five- 
dimensional cosmological constant A 5 : 


ds 2 = — u{r)dt 2 + 


2 ■ ' ,(r) 


u{r) 

u(r) = 1 — 2 M/r, 

r — 2M \ 
' - (3M/2)y 
2 M\ 


F(r) 


— Ar 


dn 2 , 


;(r) = 1 + 


F(r) 


-A 5 r 


(14) 


F(r) = 1 - 1.1241 ( —J + 1.956 ( — ) - 


2 MV 


- 9.961 




r / 
2.900 


(“) 


n 


where ' = d/dr. The function F(r) describes the per¬ 
turbation caused by the bulk. The best fit for it was 
obtained in [ 2 (| . 

The field equations induced on the brane were derived 
by Sasaki, Shiromizu and Maeda (dll : 


Gfn, — A 4 gf lv + 2 + 3 S^ v Z^, (15) 

1V1 P14 Jli P15 

where A 4 is usual four-dimensional cosmological con¬ 
stant, g^ is the metric on the brane, T is the stress- 
energy tensor of the matter localized on the brane, S^ 
is the local quadratic stress-energy correction, Z^ is the 
four-dimensional projection of the five-dimensional Weyl 
tensor. Mpu is usual four-dimensional Planck mass and 
Mpi 5 is the fundamental five-dimensional Planck mass. 

The induced metric on the brane is flat, the bulk is 
an anti-de-Sitter Space-time as in the original Randall- 
Sundrum scenario [l3|, then Z MJ/ = 0 jTgj- Therefore 
the correction term due to the contribution from ACPY 
topology m that follows from m has the form 


„ e 2 121 „ 

8 = 1 — — -M 2 , 

M e 108 


( 12 ) 


f 2 M 2 1 

a + ° (r_6) ' 


(16) 
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According to © the expansion term of the PPN-order 
should be proportional to r~ 2 . The correction (fl6l) con¬ 
tains the next perturbation order which lies beyond PPN 
similarly to the Gauss-Bonnet case ©. Therefore the 
obtained contribution (flbl) cannot be observed in the so¬ 
lar system experiments as well. This conclusion on the 
Randall-Sundrum model predictions confirms the result 
for the Figueras-Wiseman and coincides with the Gauss- 
Bonnet case. 


III. THERMODYNAMICS AND PBHS 


for GRB progenitors and therefore the limit estimation 
for black hole evaporation rate can be obtained in such 
a way. 

Different theories of gravity predict different black hole 
evaporation rate and therefore different initial masses of 
the PBHs that fully evaporate for the Universe lifetime. 
In this paper we compare the evaporation rates for the 
Gauss-Bonnet and the RSII black hole solutions. Accord¬ 
ing to the GRB data and the precision of the Fermi LAT 
telescope the closest distance d at which the telescope 
will be able to detect the evaporation of primordial black 
holes is [11| 


It is conjectured that density fluctuations in the early 
Universe could have created black holes with arbitrarily 
small masses even to the Planck scale [H}. These black 
holes are referred to as primordial black holes (PBHs) 
[36} and can be used to consider viable theories in cos¬ 
mological conditions. 

Hawking evaporation [53, HU is one of the most signif¬ 
icant properties of a black hole and can be described by 
the mass-loss rate equation (39} : 


dM 1 k B 

dt 256 7r 3 M 2 ’ 


(17) 


where M is the mass of a black hole and k B is the Stefan- 
Boltzmann constant. Hawking evaporation is a quantum 
process forbidden in classical physic. An outgoing radia¬ 
tion has to cross a potential barrier of black hole horizon 
[ 40 } so a radiation surrounding a black hole is in a ther¬ 
mal equilibrium and can be described as a black body 
radiation. Therefore the black hole evaporation obeys 
the following law: 


dM 

dt 


k B ST 4 , 


(18) 


where S is its surface area. We use this formula to esti¬ 
mate the lifetime of the black holes in the Gauss-Bonnet 
and the Randall-Sundrum models. 

According to (1171) the black holes with stellar masses 
evaporate very slowly and do not lose mass through this 
process noticeably. On the other hand PBHs with initial 
masses smaller than 


M 0 ~ 5.0 x 10 14 g (19) 

have already evaporated and can contribute to the ex- 
tragalactic background radiation [38}. PBHs with ini¬ 
tial mass greater than M 0 (fl9l) should be evaporating till 
now [4l} . According to some models of black hole evap¬ 
oration mm the last stages of this process can be 
accompanied by bursts of high-energy particles [S] includ¬ 
ing gamma radiation with energy of the MeV-TeV range 
that occur at the distances about z < 9.4 [42}. Such 
events should be rather rare and on the other hand the 
set of more simple explanations for most part of gamma 
ray bursts (GRB) exists. Nevertheless PBHs at the last 
stage of evaporation can serve as additional candidates 


d ~ 0.04 





( 20 ) 


where Q is the angular resolution of the telescope, E - 
energy range of the telescope, T - temperature of the 
black hole. Spending the same reverse procedure and us¬ 
ing a telescope detected gamma-ray bursts lead to the 
observable difference of the PBH initial mass on its fi¬ 
nal evaporation stage can vary from the GR predictions 
within the following limits: 


^investigated theory lri5 

M^r > 1 j 

Therefore we use this limit as the mass cutoff threshold 
in our calculations. 

Using the method by Shankaranarayanan, Padmanab- 
han and Srinivasan [43} it is possible to rewrite the ex¬ 
pression for the Gauss-Bonnet black hole temperature 
and then use 0. In the astrophysical case the dilatonic 
charge D cs 1/M 0- Therefore the right part of (fTSl) 
can be expanded in series as 


dM 1 k B 1 k B / ln \ 

---1---h O (M~ 10 ) 

dt ~ 256 7t 3 M 2 512 7 r 3 M 6 1 ’ 


( 22 ) 


The initial mass of the PBH that fully evaporates during 
the lifetime of the universe in this case is 


M gb = 8 x 10 14 g. (23) 

The difference between the obtained value and the similar 
GR quantity (flT)l) is smaller than the cutoff threshold 
set by (ED- Thus the specific features of the Gauss- 
Bonnet evaporation rate are negligible at the current level 
of accuracy and the predictions of Gauss-Bonnet gravity 
for the Hawking evaporation are indistinguishable from 
those of the GR. 

One of the first and most studied black hole solu¬ 
tions of Randall-Sundrum model was found by Dadhich, 
Maartens et al. [3 [3 • They obtained an exact local¬ 
ized black hole solution, which remarkably has the math¬ 
ematical form of the Reissner-Nordstrom solution, but 
without electric charge being present 0: 

2M ( q \ 1 

~ 9tt-9rr-l~— +[mTJ^2> ( 24 ) 
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The Reissner-Nordstrom-type correction to the 
Schwarzschild potential in EH) can be thought of 
as a dimensionless “tidal charge” parameter g, arising 
from the projection onto the brane of free gravitational 
field effects in the bulk transmitted via the bulk Weyl 
tensor 0. The projected Weyl tensor, transmitting the 
tidal charge stresses from the bulk to the brane, is [14(: 

= — (^ 7^/2 ^ ^4 ~ 2 TfjXv + hfi I /) , 

where = g^ + u^u^ projects orthogonal to 4-velocity 
field and r^ is a unit radial vector. 

The mass loss rate obtained similarly to the Gauss- 
Bonnet case equals to 


dM 

dt 


1 liB 
216 7r 3 M 2 


+ O (M~ e ) . 


(25) 


The leading term in (l25l) cannot produce the needed 5- 
order difference defined by the threshold parameter HUD- 
The initial mass of the Dadhich-Rezania black hole that 
evaporates completely during the lifetime of the Universe 
proves this fact: 


M D r = 5.3 x 10 14 g. 


(26) 


As the obtained difference is much less than the cutoff 
threshold m the “tidal charge” influence is vanishing 
and cannot have experimentally verifiable consequences. 

The black hole evaporation for the ACPY solution dis¬ 
cussed in the previous section also was considered in a 
similar manner. The evaporation rate of this solution 
has completely the same form as the original Hawking 
formula El up to M 10 terms thus the value of the 
initial mass equals to that given by the GR: 

Map = 5.0 x 10 14 g. (27) 


The results for Figueras-Wiseman solution are the same 
because of the form of the solution m- 

The obtained results (E2l [251 E51E71) lead to a conclu¬ 
sion that the precision of the currently existing GRB data 
is not sufficient to distinguish the GR, the Gauss-Bonnet 
and the Randall-Sundrum gravity from each other via 
the PBH consideration. 


IV. DISCUSSION AND CONCLUSIONS 

In this paper we discussed the possibilities to test the 
theories extending GR in different ways by the example 
of the Gauss-Bonnet and the Randall-Sundrum models 
both in the weak field and the cosmological limits. For 
this purpose the post-Newtonian expansion and the black 
hole evaporation of these theories’ solutions were consid¬ 
ered. 

As the Gauss-Bonnet term coupled with the scalar field 
does not influence the post-Newtonian limit © there¬ 
fore the nontrivial scalar hair generated by it [7], 11] does 


not contribute into the required order of the spherically- 
symmetric solution’s expansion. This result agrees with 
the previous conclusions by Sotiriou and Barausse [29j 
who considered the cosmological solution of the action 
© and showed that the influence of the Gauss-Bonnet 
term is negligible at solar system scales. Combining these 
two results we can state that the leading term of Love¬ 
lock expansion © describing second-order curvature cor¬ 
rection does not provide any visible deviation from GR 
predictions in the weak-field limit therefore such theory 
of gravity fully agrees with GR. 

This conclusion is also valid for any model with higher- 
order curvature corrections having a proper Newtonian 
limit. As the Gauss-Bonnet term is the leading curva¬ 
ture correction of the Lovelock gravity its contribution 
into the post-Newtonian expansion of the metric also is 
the largest one. Taking into account the others Euler 
characteristics i.e. the next orders of curvature correc¬ 
tions isn’t able to change the picture as their influence 
is even less and obviously lies far beyond the PPN limit. 
Thus the conclusions for the Gauss-Bonnet model can be 
generalized to the Lovelock gravity. 

The theories with curvature power series are not the 
only method of geometrical extending of GR. In the 
generic case the Lanrangian can contain an arbitrary 
function of the Ricci scalar R. Such theories set up 
f(R )-gravity [3 |45| and Lovelock gravity is its partic¬ 
ular case. Many /(/?)-gravity models such as ln(R) or 
1/R [3 13] were originally introduced as attempts to 
explain dark energy or dark matter. They don’t have a 
proper PPN limit 44] and are inapplicable to solar sys¬ 
tem scale. Therefore our conclusions for Gauss-Bonnet 
theory in the weak-field limit are applicable for Lovelock 
gravity and /(R)-gravity of the Lovelock type. 

The thermodynamical properties of the Gauss-Bonnet 
black hole solution were considered in details earlier 
however only the black holes of Planck scales 
were investigated. For the black holes with larger mass 
the influence of the Gauss-Bonnet term and the scalar 
field becomes negligibly small therefore the evaporation 
is predictably the same as in the GR case. 

Since Randall and Sundrum proposed a theory of grav¬ 
ity with noncompact extra dimension 1/ 
black hole solutions have been found [U 
Consideration of the post-Newtonian expansion of the 
Figueras-Wiseman solution [3 reveals such possible ef¬ 
fect as a negative nonlinearity of gravitational superpo¬ 
sition (fI51) . It naturally results from the theory itself 
because gravity is allowed to propagate to the extra di¬ 
mension in Randall-Sundrum model. However the break¬ 
ing of gravitational superposition turns out to depend on 
a negligibly small parameter thus the predictions of the 
Figueras-Wiseman solution fully agrees with GR and the 
present observations. This effect may influence the strong 
field regime (close binary systems, black holes) as a con¬ 
sequence of curvature growth. So the next step could be 
the search of such features of the Randall-Sundrum model 
in the strong field limit. Fortunately this investigation 
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is admissible as the large stable black hole solutions for 
RSII black holes have been found [H, [26|. 

The consideration of the other black hole solution by 
Abdolrahimi, Page et al. [26j| shows that the terms de¬ 
scribing the bulk influence (fThl) greatly exceed the limits 
of the post-Newtonian approximation. As a result both 
large Randall-Sundrum black holes solutions cannot be 
distinguished from the Schwarzschild metric at the solar 
system scales. 

We have examined the evaporation rate for the 
Randall-Sundrum black holes as well. The results for 
one of the first solutions obtained by Dadhich, Maartens, 
Papadopoulos and Rezania [l4] and the latest one by Ab¬ 
dolrahimi, Page et al. [26j are presented (|22l l23l I25ll27l . 
The difference between the Dadhich-Rezania solution and 
the GR is negligibly small and the Page solution coincides 
with GR completely. 

As is easy to see, many extended gravity models cannot 
be distinguished from GR and from each other neither at 
the solar system scales nor by the black holes thermo¬ 
dynamic properties. Therefore the coincidence of these 
extended theories with the GR serves a good argument 
in favor of their validity. However it does not mean that 
no difference can be found by other verification meth¬ 
ods. Except the weak field and the cosmological tests a 
strong field approximation is widely used. It has such a 
verification laboratory as close binary systems first of all 
containing pulsars as one or even both its components. 
A great amount of data has been obtained from these ob¬ 
servations and it obviously should be used for testing the 
extended gravity models though this method has its own 
shortcomings. If the orders of the post-Newtonian correc¬ 
tions of the extended gravity models lie beyond the PPN 
order it is natural to suggest that the parameters of the 


models should be limited via the second or the third post- 
Newtonian orders. The corresponding formalisms of the 
2PN and 3PN really do exist [13,S3- These formalisms 
consider the gravitational radiation and its subtle effects 
on pulsar timing and orbit parameters. However many 
calculations there are based on GR and do not suit for 
comparing arbitrary extended theories of gravity like the 
PPN formalism [27] does. 

There are also the other ways to test astrophysical 
predictions of the extended theories of gravity such as 
accretion onto massive objects and microlensing. After 
computing the accretion rate for some solution the result 
can be compared with GR predictions and some other 
extended gravity cases. The investigation of the data 
of gravity lensing events also is a perspective method as 
these data become more and more complete. Verifica¬ 
tion of the extended gravity models via studying binary 
systems and particularly the pulsar data needs special 
methods and approaches. Their construction is the sub¬ 
ject of further considerations. 
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